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SMOOTH BOUNCE IN AFFINE QUANTIZATION OF BIANCHI I 


HERVE BERGERON, ANDREA DAPOR, JEAN PIERRE GAZEAU AND PRZEMYSLAW 

MALKIEWICZ 

Abstract. We present the affine coherent state quantization of the Bianchi I 
model. As in our previous paper on quantum theory of Friedmann models, we 
employ a variable associated with a perfect fluid to play a role of clock. Then 
we deparameterize the model. A distinctive feature, absent in isotropic models, 
is an extra non-holonomic constraint, which survives the deparameterization and 
constrains the range of physical variables. The appearance of the constraint reflects 
the ‘amplification’ of singularity due to anisotropy. The quantization smoothes 
the extra constraint and allows quantum contracting trajectories to be smoothly 
transformed into expanding ones. Making use of affine coherent state we develop 
a semiclassical description. Figures are included to illustrate our result. 


Contents 

1. Introduction 2 

2. Spacetime geometry 3 

3. Canonical framework 4 

4. Quantization of true Hamiltonian h 6 

4.1. Affine quantization: a compendium 6 

4.2. Quantization of Hamiltonian 7 

5. Analysis of positivity constraint h > 0 8 

5.1. Positivity via operator modification 9 

5.2. Positivity via domain restriction 9 

5.3. Positivity via redefinition of classical variables 9 

5.4. Positivity via redefinition of classical expression 10 

6. Refined quantization of true Hamiltonian 10 

6.1. General formula for the operator 10 

7. Semiclassical approach 11 

7.1. From Lagrangian 11 


Date : February 2, 2015. 

P.M. was supported by MNISW Grant “Mobilnosc Plus”. A.D. was supported by the grant of 
Polish Narodowe Centrum Nauki nr 2013/09/N/ST2/04312. The work was supported by Polonium 
program No. 27690XK Gravity and Quantum Cosmology. 

1 



2 


H BERGERON, A DAPOR, J-P GAZEAU AND P MALKIEWICZ 


7.2. Lower symbol 

8. Singularity resolution 

8.1. Numerical examples 

8.2. Lower symbol of the Hamiltonian 

8.3. Effective dynamics of q and p 

8.4. Effective spacetime 

8.5. Effective dynamics of volume and expansion rate 

9. Discussion 

9.1. Mechanism behind the singularity resolution 

9.2. Features of the quantum bounce 

Appendix A. Computation of quantized Hamiltonian 

Appendix B. Quantization of 9 

References 


12 

12 

12 

13 

14 

15 

15 

16 
16 

17 

18 

19 

20 


1. Introduction 


In our previous paper [1] we proposed a quantization of spatially homogeneous 
and isotropic spacetime, based on coherent states for the affine group (also known 
as wavelets). The main result of this approach was the appearance of a quantum 
centrifugal potential, whose effect is the regularization of the singularity, replacing 
the Big Bang with a bounce both at quantum and semi classical levels. Moreover, 
the true Hamiltonian is essentially self-adjoint, and hence the dynamics is unitary 
and unambiguous. The results obtained in this simple model of cosmological collapse 
encouraged us to tackle further investigations concerning more elaborate models of 
the Universe. 

In the present work, we study the anisotropic case of Bianchi I, allowing for more 
degrees of freedom in the description of spacetime. The canonical classical variables 
p and q introduced in [1] apply equally well to the present case, so that we can make 
contact with our previous result. It needs to be said that we are entering an area 
much less visited than the Friedmann models, and thus we should expect some new 
difficulties. 

In the FRW case [1], apart from studying the quantum true Hamiltonian, we were 
also interested in the properties of quantized volume and expansion rate, as they 
signal the singularity by vanishing and blowing up respectively. While the quantum 
volume turned out to be self-adjoint, the quantum expansion rate did not. Never¬ 
theless, it was possible to give it a semiclassical meaning through its lower symbol 
(i.e., its expectation value on affine coherent states, opportunely peaked at classical 
phase space points). In the present work, in addition to these quantities, we also 
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have to deal with shear and distortion describing the anisotropic evolution. Classi¬ 
cally, the shear and some of the scale factors blow up at the singularity. However, 
the existence of minimal volume for the universe at the quantum level prevents the 
unbound growth of anisotropy, i.e. it smoothly evolves through the bounce. We will 
not provide an explicit evolution of shear and distortion at the quantum level, but 
confine ourselves to giving a well-behaved semiclassical Hamiltonian, which generates 
a well-behaved dynamics of anisotropic variables. 

We should point out that, as in any Hamiltonian quantum theory of cosmology, 
we must choose a degree of freedom to play the role of physical time. (In the reduced 
phase space scheme, this is done before quantization.) Following Schutz [2, 3], we 
identify this internal clock with a barotropic fluid filling the Universe. This partic¬ 
ular choice is technically attractive because the inclusion of fluid extends the phase 
space with canonical variables T and px, where the fluid’s momentum px is a con¬ 
stant of motion. Thus, it is possible to perform a reduction to the physical phase 
space equipped with a time-independent true Hamiltonian. This feature may be 
particularly useful in future analysis of non-integrable Bianchi IX model. 

In Section 2 is briefly presented the geometry of the model we wish to quantize. 
Its canonical formulation is described in Section 3, where the true Hamiltonian h is 
defined, and the non-holonomic constraint h > 0 is specified. The affine coherent 
state quantization of the true Hamiltonian is performed in Section 4 without tak¬ 
ing into account this constraint. In Section 5 we discuss different possibilities to 
implement the latter. It appears that the best approach consists in incorporating 
this constraint in the classical true Hamiltonian to be quantized. To do this, we use 
the procedure for coherent state quantization of distributions presented in [4], The 
resulting quantum Hamiltonian is presented in Section 6 (detailed calculations can 
be found in Appendix A), and its lower symbol is computed in Section 7. This can be 
thought of as an effective classical Hamiltonian accounting for quantum corrections. 
An approximate analytic expression for this effective Hamiltonian is used in Section 
8 to study the effective dynamics of fundamental variables, volume and expansion 
rate. Numerical analysis of the Big Bang singularity resolution is also presented. 
Section 9 concludes with a discussion of some interesting issues of our results and 
possible future directions. 


2. Spacetime geometry 

In the present work we study the anisotropic evolution of compact (say, with torus 
topology T 3 ) locally flat spatial sections of the spacetime equipped with the metric 

ds 2 = — N(t) 2 dt 2 + ^2 di(t) 2 (dx 1 ) 2 . 


( 1 ) 
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We fix the spatial coordinates via J gl dx l = 1. The Ricci scalar reads: 

4 -a ^\n) ~\mrj_ ■ 

Filled with perfect fluids, these spacetimes are singular as the world-lines of co¬ 
moving observers terminate within a finite proper time. Starting from a spherical 
chunk of space, the anisotropic evolution will deform its initial shape. Following this 
observation, the spacetimes are said to have cigar or barrel type singularities, with 
two scale factors vanishing and the other one blowing up or two scale factors reaching 
a finite value and the other one vanishing, respectively, at the singularity. A more 
detail discussion of solutions may be found in [5]. 


(2) R = 2 


3. Canonical framework 

The Hamiltonian constraint of general relativity in the ADM variables (<?p,p*' ? ), 
equipped with the Poisson bracket { qij(x),p kl (x ')} = 167T GS^S^8(x — x'), reads [6] 


( 3 ) 


where 


(4) C° 


H, ; = ^J [NC° + N t C n ] d 3 x , 


-\/det(ftj) 


LR + 


det(gjj) 


2 (Pk) 2 -P ij Pa 


, C' = -2p% 


are first-class constraints and N, Ni are Lagrange multipliers and k = 167tG. The 
truncation of the gravitational variables to the cosmological sector translates into 
qij(x) := 5ijO% and p^{x) := 8ijp\ which in turn lead to C l = 0 and 


( 5 ) 


H, = — 


1 N 


K, d\Qj2&Z \ 2 


<S ^ ( hV V ? ~ 


2i\2 


where we defined 

i 1 (02^3) 2 1 (« 3 a i) 3 1 (01O2) 

6 P = ~Tr -> P = ~m -> P = — Tt -> 

N ai N a 2 N a 3 

satisfying {a?,jp} = nSP . The Hamiltonian (5) can be brought to a simpler, diagonal 
form. We achieve this by Misner (canonical) parametrization of the phase space [7]: 
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canonical pairs are given by (/3°,po, /3 + ,p + , (3 ,p~), where 


( 7 ) 

( 8 ) 



/ 1/3 
1/6 
1 

\2\/3 

a 


1/3 1/3 \ 

1/6 -1/3 

1 

2\/3 ° / 

1 1\ (2a\ Pl \ 

1 -2 2a|p 2 

-x/3 0 / \2a 2 3 p 3 J 


( lnaA 
lna 2 
In a 3 ) 


? 


In term of these variables, the gravitational Hamiltonian reads 

( 9 ) = (-Po+pI+P-) ■ 

At this point, H g is still a constraint. To make it into the physical true Hamiltonian, 
we include matter in the form of a barotropic fluid subject to the equation of state p = 
wp (w = const). Following Schutz [2, 3] we obtain for this matter the Hamiltonian 

( 10 ) H m = N-^ = Ne-^»°p T , 

where T and pt are canonical variables associated with the fluid. In analogy with 
our previous work on the FRW models [1], we now replace (/3°,po) with the following 
canonical pair: 


( 11 ) 


3(1—uj) a 0 

q := e 2 v , p ;= 


3(1 — w) 


_ 30-^0 


Po 


such that [q\ = 1, [p] = Lr and {q,p} = k. At the physical level we notice that this is 
exactly the same pair of observables which we considered for the FRW models, i.e., 


( 12 ) 


q = (0102^3) 2 ™, p = 


3(1 


w 


-(dia 2 a 3 ) 2 0, 


where 0 is the expansion rate. The full Hamiltonian constraint reads 


(13) 



0, 


where p := —2w/(l — w), a := 3(1 — w)j2 > 0. Solving the constraint for px leads 
to 


( 14 ) 


«pt = ^p 2 - ^(pI + p-)q 2 
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Upon choosing T/k as the clock variable, we obtain the true (non-vanishing) Hamil¬ 
tonian 

(15) h = ^p 2 - ^(p+ +p 2 _)q~ 2 . 

describing the evolution of gravitational degrees of freedom (q,p, /3~,p-, /3 + ,p + ) G 
x M 5 . The most important difference between (15) and the true Hamiltonian 
found in [1] is the appearance of the attractive potential proportional to —q~ 2 - We 
must then impose positivity of h as a separate, non-holonomic constraint 

(16) h > 0 , 

which classically restricts the range of gravitational variables. As a side remark, 
we notice that in the limit of vanishing anisotropy, p- —> 0 and p + —» 0, the true 
Hamiltonian of the flat FRW spacetime [1] is recovered. 

4. Quantization of true Hamiltonian h 

4.1. Affine quantization: a compendium. As the complex plane is viewed as 
the phase space for the motion of a particle on the line, the half-plane is viewed as 
the phase space for the motion of a particle on the half-line. Let the upper half-plane 
H+ := {(g,p) | p G M, g > 0} be equipped with the measure dqdp. Together with 
the multiplication 

(17) {q,p)(qo,Po) = {qqo,po/q + p), q g R+, p g m, 

the unity (1,0) and the inverse 

(18) (g,p) _1 = , 

n + is viewed as the affine group Aff + (M) of the real line, and the measure dqdp is 
left-invariant with respect to this action. The affine group Aff + (M) has two non¬ 
equivalent unitary irreducible representations (U1R). Both are square integrable and 
this is the rationale behind affine integral quantization [4], The U1R U + = U is 
realized in the Hilbert space V. = L 2 (M^_, dx): 

(19) U (q, p)ip(x) = (e ipx /y/q)ij(x/q). 

Picking an admissible unit-norm state ip, i.e., an element in G L 2 (M^, da;)nL 2 (M^, dx/x), 
named fiducial vector, produces all affine coherent states (or wavelets) defined as 

( 20 ) | q,p) — U(q,p)\ip). 

These states in (q, p)-representation, i.e., as functions ( q ,p(q',p') ■= {q',p'\q,p), are 
expected to be well peaked around their phase-space point (q,p). 
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Square integrability of the UIR U and admissibility yield the resolution of the 
identity 

/ \ f dqdp dx o 

( 21 ) ^r~Jq,p)(q,p\ = I, where c 7 := J . 

The covariant integral quantization of classical observables f(q,p) follows: 

(22) / A f = [ f(q,p)\q,p){q,p\p^~ ■ 

Ju+ Z7TC_i 

The above map is linear, transforms the constant function 1 into the identity oper¬ 

ator, and real functions into symmetric operators (provided that the integral makes 
sense). Also, note that this quantization procedure can be easily extended to singular 
functions and distributions [4], In the sequel we restrict the choice of fiducial vector 
to real functions. 

The map (22) yields canonical commutation rule, up to a scaling factor, for q and 
p. Indeed, in representation “ x ”, 

(23) A p = = P,A q = (co/c_i) Q , Qf(x ) := xf ( x ) , [A q , A p ] = (co/c_i)iJ. 

Multiplication operator Q is (essentially) self-adjoint whereas P is symmetric but 
has no self-adjoint extension [ 8 ]. The property of the procedure which is crucial for 
our present purpose is the regularization of the quantum kinetic energy: 

(24) A p 2 = P 2 + q^Q~ 2 with ^ := f — u (u)) 2 . 

Jo c -i 

The additional term in regard with the standard canonical quantization is a cen¬ 
trifugal potential cx Q~ 2 whose strength depends on the fiducial vector only, and can 
be made as small as one wishes through an appropriate choice of ip. If we consider 
the quantum dynamics of a free motion on the open half-line M+, it is proved [ 8 ] 
that the operator P 2 = —d 2 /dx 2 alone in L 2 (M^_,da;) is not essentially self-adjoint 
whereas the regularized operator (24) is for ^ > 3/4. It follows that for < 7 ^, > 3/4 
the quantum dynamics is unitary during the entire evolution, in particular in the 
passage from the motion towards the origin (contracting branch) to the motion away 
from it (expanding branch). 

4.2. Quantization of Hamiltonian. We quantize the (q,p) variables introduced 
in ( 11 ) according to the above affine quantization framework, while for the other 
variables we choose the canonical quantization one since they respective phase space 
geometry is the harmless plane M 2 . Making use of the formulas above (see also [1]), 
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the quantum version of the true Hamiltonian (15) is 


(25) 


h >-4 ylh — —P 2 + 2 - + P~)Q 2 


The operators P + and P_ are the usual momentum operators on the real line acting 
in L 2 (M 2 , d/3 + d/3 — ). Since they are self-adjoint and commute with the Hamiltonian, 
we can replace them with some fixed real eigenvalues, k + and fc_, which will remain 
constant during the evolution. 

As in [1], we choose ip as the function smooth in (0, oo) and rapidly decreasing at 
0+ and at the infinity: 

(26) ip = ip u (x) = — e ~ 3 (^ 12+ ^n), with v > 0 and £ [2 = > 0- 

y/2xK 0 (v) A(ffi) 

Here and in the following, K r {z) denotes the modified Bessel functions [9]. Since we 
deal with ratios of such functions throughout the sequel, we adopt the convenient 
notation 


(27) 


P,rs P,rsiy) 


K r (u) 

K a {u) 


1 


Ssr 


One convenient feature of such a notation is that £ rs (i/) ~ 1 as u -A oo (a consequence 
of K r {v) ~ \J 7r/(2z/)). Hence, the quantum Hamiltonian (25) now reads 


(28) 


A = ^P 2 + a 2 P ^(v)Q 2 - (k% + k 2 _)Q 2 


where s(v) = q^v = (1 + z/^oi) /4. It is convenient to introduce the auxiliary strength 
parameter 

(29) := q(v) - 

a z a P 

so (28) reduces to 

a 1 

(30) A = — (P 2 + x{v) a 2 p Q 2 ) • 


The key difficulty is the additional non-holonomic positivity constraint (also present 
in open [k = —1) FRW models), h > 0, to which the next Section is dedicated. 


5. Analysis of positivity constraint h > 0 

In what follows we consider different possible ways of implementing the positivity 
constraint (16) in the quantum theory. 
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5.1. Positivity via operator modification. The positivity requirement of the true 
Hamiltonian may be imposed directly on the operator A h by just putting x(z/) > 0. 
This leads to the condition on the amount of shear, which is expressed in terms of 
eigenvalues k + and 

(31) k\ + k 2 _ < a 2 a 2 P (£ i2 ) 2 £i-i<T- 

The above condition reduces the space of allowed quantum states by restricting the 
spectral values associated with the shear operator. More precisely, the amount of 
anisotropy is bounded from above by the right hand side, meaning that a highly 
anisotropic spacetime cannot be realized. However, there is no reason why the Uni¬ 
verse should be so much isotropic. On the contrary, since classically the anisotropy 
goes as ~ a -6 with the overall scale factor, one reasonably expects that the Universe 
is more and more anisotropy-dominated the closer it is to the classical singularity. 
We thus conclude that the above condition for the allowed amount of anisotropy in 
the Universe is too restrictive, and discard this approach as unphysical. 


5.2. Positivity via domain restriction. Let us consider the eigenvalue problem 
for operator A h in (30) on L 2 (R* +1 dx): 


(32) 


d 2 i/> x{v) a 


x z 


0 = A 0 


A possible way of implementing the condition (16) is to ensure the self-adjointness 
of (32) and then constrain the Hilbert space to the positive eigenvalues, A > 0. The 
difficulty in this approach is to implement such a Hilbertian projection in analytical 
computations, since it requires full knowledge of the spectral decomposition of Ah. 
We drop this approach for technical reasons. 


5.3. Positivity via redefinition of classical variables. Consider the following 
canonical transformation: 


(33) 

(34) 


, 1 

q i—^ q - 


pq 


2 p 2 — p 2 /q 2 ’ 


/ 2 b 

p^p -=p --, 

qz 


+ In 

2a-p \pq + p 


p±^p ± :=p± 


where p 2 : = (p++p 2 ) /° 2 - In these new variables, the positivity constraint (16) 
becomes p' > 0. There is also another constraint, which is \q'p'\ > p/2. We observe 
that this condition splits the physical phase space into two disconnected regions. 
If we now define q' := q'p' — p/2 and p' := hip', then ( q',p ') £ x R, and 
the true Hamiltonian reads h = . The advantage here is that we do not have 

to impose the positivity constraint quantum mechanically, since h in this form is 
automatically positive-definite. The disadvantage of this approach is that we lose 
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connection with the variables used for FRW models [1], More importantly, the 
expanding and contracting branches of the Universe become disconnected. For these 
reasons we discard also this approach. 

5.4. Positivity via redefinition of classical expression. Another idea is to take 
0(h)h (where 9 is the Heaviside function) as the classical true Hamiltonian, and 
quantize it in the spirit of integral quantization framework [4, 10]. The resulting 
operator Aq^ will have the positivity constraint implemented. Moreover, in this 
way we do not redefine the initial Hilbert space. In the next Section we follow this 
last idea, in what we call “refined quantization” (as opposed to Section 4). 


6. Refined quantization of true Hamiltonian 


6.1. General formula for the operator. Let us consider the simplified form of 
the Hamiltonian (15) 

(35) h s = p 2 - 

q 

obtained after dropping unnecessary factors, together with the positivity constraint 
h s > 0. We derive from the decomposition h s = 0(h s )h s + 9(— h s )h s and the linearity 
of the affine integral quantization the equation 


(36) A 0 ( hs)hs - A hs - A 0( _ hs ) hs . 

With the fiducial vector (26) and from results given in [1], we obtain for the affine 
quantized version of h s in u x” representation: 


(37) 


(UA.M = 


-^2 + (^)- 6 l O ') 2 * 2 ) ^2 


5(x — x') 


All the details about the derivation of this expression are given in Appendix A. 

The second term, Ae(_h s )h s for the same fiducial vector (26), in “x” representation 
reads 


(38) 


(z'l-V-MhJz) 


2 

7riF 0 (z/)c_i 


( fcRe[^ 0 (2 7 )] 

\\Zxx'(x — x') 2 


4£i2 Vxx i lm['yK 1 (2'y)\ \ 
v (x — x') 3 (x + x') ) ‘ 


where 


(39) 



4 k“ 

i - —(x — x')) 

Ul2 



4 k 

Ul2 


( 1 - 2 /)) 


with y = x'/x. This is a non-local operator, and there is little chance to solve its 
eigenvalue problem analytically. We are not going to analyze it except for noticing 
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that the operator A 0 ( hs ) hs is positive-definite and as such it admits self-adjoint ex¬ 
tension^), e.g., Friedrich’s extension. Therefore, the existence of unitary evolution 
is guaranteed. In what follows we make use of the semiclassical description available 
in our approach based on affine coherent states. 

7. Semiclassical approach 

7.1. From Lagrangian. Inspired by Klauder’s approach [11, 12], we present a con¬ 
sistent framework allowing to approximate the quantum Hamiltonian and its associ¬ 
ated dynamics by making use of a semiclassical Lagrangian approach. 

For a general Hamiltonian operator H, the Schrodinger equation, ih^\^!(t)) = 
hi m), can be deduced from the Lagrangian: 

(40) L(T, T,A0 := <tt(t)| (ih^ - h) | *(t)), 

through the application of the variational principle with respect to |\k(f)). Following 
Klauder, we assume that |(t)) is an affine coherent state \q(t),p(t)), where q(t) 
and p(t) are some time-dependent functions. Then the Lagrangian (40) turns to the 
semiclassical form 

L (q,q,p) = (q(t),p(t )| \q(t),p(t)) 

(41) = -£02 {y)qp ~ ( q(t),p(t)\R\q(t),p(t )) . 

We rescale the family of coherent states to define a new family of coherent states 
| q(t),p(t)) new := |A q(t),p(t)), where A = I 11 this way, we ensure that (q,p), 

computed as the expectation values wrt to the rescaled family, fully correspond to 
the classical pair, i.e. (q,p) = (q,p)- Then, we derive the semiclassical equations of 
motion as 

(42) q ( ^Hu/.//). 

(43) p = —^H(q,p). 

where H (q,p) is the lower symbol of H. H (q,p) = (\q(t),p(t)\tl\\q(t),p(t)). We note 
that the classical relation between q and p cannot hold any longer on the semiclassical 
level since this relation is now given by (42), in which classical H is replaced with 
H. This is viewed as the consequence of the quantum non-commutativity of basic 
variables. In our reconstruction of the semiclassical description of spacetime we will 
keep the interpretation of q as fundamentally given by (12). 

We now specify H = H 0 ( hs ) hs , and H = (Ag,p|H<9( hs ) h JAg,p). The latter lower 
symbol is the basis of our semiclassical description. 
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FIGURE 1 . The figure shows the trajectories generated by the lower sym¬ 
bol of the true Hamiltonian. The trajectories are clearly reversed and the 
bounce occurs. The parameters are chosen as k = 10 ap (for the potential 
strength), v = 3 (for the fiducial vector). 


7.2. Lower symbol. We find that the lower symbol of Hg(_h s )h s reads 

2 

{Xq,p\A e{ -h s )h s \Xq,p) = - 


Re [K 0 (27)] _ 4£i 2 Im[7/ii(27)] \ 7^(27) y 2 d y 
y{i-y) 2 v (1 - y) 3 (y + 1)/ (y 2 + 1) 2 c_ x 

where 

and y = x' jx. The analytic evaluation of the above integral is rather intractable and 
we proceed with a numerical integration instead. 

8. Singularity resolution 

8.1. Numerical examples. The semiclassical refined true Hamiltonian is defined 
as the difference of the lower symbol of (37) and the lower symbol (44). Figures 1 
and 2 illustrate the singularity resolution in the half-plane ( q,p) for different values 
of the strength k. The contour plots of the lower symbol of the Hamiltonian contain 
semiclassical trajectories as const ant-value levels. The numerical computations prove 
that they exhibit bouncing behavior. 
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FIGURE 2. Here k = 50 ap, the semiclassical trajectories are initially more 
divergent, but they are again reversed. On the left, the same scale as above. 
On the right, the plot is zoomed in, so that we can see the reversal of 
trajectories near q = 0. 


8.2. Lower symbol of the Hamiltonian. Combining the lower symbol of the 
operator A^ s , given in Appendix A, with (44), we find that the lower symbol of the 
true Hamiltonian (36) has the following form: 

(45) (Xq,p\A e(ha)hs \Xq,p) = ( p 2 X 2 q 2 + A{v) - B(v)k 2 + F u (k 2 ,p 2 X 2 q 2 )) , 

where 

(46) A{v) := -p £,o £12 ( ,y £32 — 1) , B(v) := ^20 , 
and 

(47) j^F u (k 2 ,p 2 \ 2 q 2 ) := -{\q,p\A e{ _ hs)hs \\q,p} 

is the corrective term due to the positivity constraint, for which we find the following 
limits: 

(48) F u { 0, 0) = 0 , lim F u (k 2 , 0) = — A(u) + B(u)k 2 , lim F u (k 2 , x) — 0 . 

fc—>-oo X^fOO 

Note that the lower symbol (45) is even in k and in pq. 

Based on our numerical simulations, we can guess the following approximative 
form for F u (k 2 ,p 2 X 2 q 2 ): 


( 49 ) 


F„(k 2 ,p 2 \ 2 q 2 ) ~ a u (k 2 ) b„(k 2 ,p 2 \ V) 
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FIGURE 3. (A) Plot of the exact F^(k 2 ,0) in blue and its approximation 
a^{k 2 ) in red. (B) Plot of the exact F p^ k f^ in blue and its approximation 

b^(k 2 ,x 2 ) in red. Plotted for k = 1,2,3,4. Beyond the plotted range of 
x the effect of theta function becomes negligible due to the relatively large 
positive value of {Xq,p\A e ^^ s \Xq,p), and the accuracy of the approximation 
is not very significant for reproducing the correct dynamics. 


where 

(50) 

(51) 


a u (k 2 ) = -A{v) 
b u (k 2 ,p 2 X 2 q 2 ) = 


X\{u)k 2 


+ B{v)k 2 , 


1 + Ai (u)k 2 

1 + Ai {v)k 2 

1 + Ai (u)k 2 + A 2 (z/)p 2 A 2 g 2 


That this factorization is indeed a good approximation of F u (k 2 ,p 2 X 2 q 2 ) is confirmed 
for v = 3 by the plots in Figures 3a and 3b with Ai ~ 0.3 and A 2 ~ 0.7. 


8.3. Effective dynamics of q and p. Let us now restore the original physical 
factors a and ap present in (25). The approximate form of the lower symbol (45) of 
the true Hamiltonian now reads 

( 52 ) 

(Xq,p\A g(h)h \Xq,p) « [p 2 X 2 q 2 + A{v) - B{u)^ + a v b v (^,p 2 X 2 q 2 ^ 

= h (q,p) ■ 

Making use of the approximation we obtain 

(53) (pqj = {pq, h} = 2h > 0 
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where h, our semiclassical Hamiltonian, is constant and thus 
(54) pq = 2h (T - T bounce ). 

We set T boU n C e = 0. Now we easily integrate (42-43) and obtain for q the expression 


(55) 


9 = 


' q 2 h 


1 a 2 


h 24A 2 


^4h 2 T 2 A 2 + A{v) - B(u)^ + a u (^j b 


1/2 

— + ( — I I -,4h=T^2 


>Yh 


or 




where we substituted pq according to (54) and where h on the right-hand side is 
treated as a constant of integration. Next, we find for p 


(56) p = 


2h T 


■ S’ (4W 2 A 2 + AM - B( v )£ + a„ (i|) K ( g,4h^ A 2 


k 2 


8.4. Effective spacetime. From the initial definitions given in Section 3 we have 

(57) ds 2 = — iV 2 dT 2 + ( aidx l ) 2 + ( a 2 dx 2 ) 2 + ( a 3 dx 3 ) 2 , 
where 

(58) N = — —, ai = g«e /3++ ^ , a 2 = q^e^ + ~^ , a 3 = q^e~ 2,3+ . 

qv 

In order to construct the semiclassical spacetime, in addition to q(T), we need the 
semiclassical dynamics of /3 , which can be derived from the respective Hamilton 
equations after replacing k 2 with p 2 + + p 2 _ in (52). However, the resulting formulas 
are involved and will be omitted. Instead we focus on the overall expansion and the 
total volume. 


8.5. Effective dynamics of volume and expansion rate. The volume V and 

expansion rate 0 read respectively: 


(59) 


V(T) = q- 




Towards the singularity, the volume decreases until it reaches its minimal value, 
which is 


(60) 


V m 


A(is ) a 2 \ 2a 

l + A 1 g24A 2 h y ) 


We note that for k 2 = 0 the above formula gives the minimal volume of the flat 
FRW universe [1]. The precise moment of the maximal value of expansion rate, Tq, 
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is very difficult to obtain. Let us assume that in the vicinity of the bounce, where the 
contraction reaches it maximal value, we may apply the following approximation: 


(61) 

Then, we find 


\ 


cr 


h 24A 2 


( 


4h 2 T 2 A 2 + 


A{v) 


1 + X\{v)K 


(62) 

and 

(63) Q max 


T e ~ 


A(y) 


4A 2 h 2 (l - At) 1 + A!(i/)S ’ 


0(T e ) « ±| 


/n — 2 a 2 A{y) \ 2 
^-!24hA 2 l + A 



A 4 1 + Ai(i/)^a 


9. Discussion 

9.1. Mechanism behind the singularity resolution. In the present paper we 
have dealt with the anisotropic cosmological model of Bianchi I. Classically, such 
spacetime presents a singularity which is much stronger than the one for the Fried¬ 
mann models considered in our previous paper [1], The classical trajectories in the 
(q,p) plane diverge as they approach the singular point. Indeed, the contracting and 
expanding branches are disconnected manifolds of the constraint surface. Therefore, 
one does not expect that canonical quantization may resolve such singularities. 

We employ a more suitable framework, known as the affine coherent state quanti¬ 
zation (ACS). In order to resolve the singularity with a bounce, we introduce a phase 
space in which the expanding and contracting branches are connected by a classi¬ 
cally forbidden region. Classically, an extra positivity constraint is present, to ban 
that region. The ACS quantization, because it can also be applied to distributions, 
smoothes this constraint, allowing the semiclassical trajectories to cross the classi¬ 
cally forbidden region and join the contraction and expansion with a smooth bounce 
(see Figure 1). Thus, we have obtained an anisotropic singularity resolution by a 
new mechanism, which is peculiar to our quantization framework. To illustrate the 
smoothing which takes place in the phase space, we analyze the case of theta function 
in Appendix B. Upon quantization, the theta function becomes a positive operator 
with a trivial kernel and its lower symbol is presented in Figure 4. Consequently, all 
quantum configurations corresponding to the phase space become accessible by the 
physical (i.e. satisfying all the classical and quantum constraints) motion. 

The present result shows that, in order to solve the singularity problem, one needs 
(at least partially) to impose the gravitational constraints at the quantum level. 
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Figure 4. Contour plot of the semiclassical expression (q,p\Ag\q,p). 
Here k = 3, q is running on the interval [0,2] and p on [—10,10]. 
(q,p\Ag\q,p) is running in the interval [0.74,1], 


It does not mean that the Dirac approach is preferred. Rather, a combination of 
reducing the constraint partially on the classical level and partially on the quantum 
level seems to be the best option. 

9.2. Features of the quantum bounce. We notice in formula (60) that the larger 
the anisotropy k 2 , the smaller the minimal volume which is reached by the collapsing 
universe. This straightforward observation challenges the relevance of the Planck 
scale to the big bounce as the change in anisotropy may make the bounce scale 
arbitrarily low. 1 Moreover, we notice that in (60) the minimal volume depends only 
on one effective constant A(u), that is, the one which was present already in the 
isotropic models. Thus, we have obtained an extension of our previous result [1], 
which is exactly the “isotropic limit” of our present result. The next step in our 
research is to study in similar manner other types of singularity. In farther future, 


1 Another argument against the role of Planck scale in dynamics of quantum cosmological models 

stems from the lack of privileged time function in quantum gravity. A detailed reasoning may be 
found in [13]. Even for a specific choice of time, it has recently been shown [14] that quantum 
gravitational effects could be detected by relatively low-energy particles, provided that the spacetime 
is in a sufficiently non-classical state. 
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we will consider the effect of the obtained singularity resolution on cosmological 
perturbations. 

Appendix A. Computation of quantized Hamiltonian 

The affine CS quantization (with fiducial vector if;) of the simplified Hamiltonian 
h A , = p 2 — k 2 /q 2 together with the positivity constraint h s > 0 is performed through 
the decomposition into two terms 

A6»(h s )h s = A hs - A0(_ hs ) hs , 

where the operators read in “x”- representation, (x'\A\x) = A, 

4 f d ?d Pnf-2 k ' 2 \ ( 2 k 2 \ !-///>../ / , 


Ah. = 


9 \P 2 - “2 ) \P 2 - — ) e%V{X X ' )( l ^{x'/q^ix/q ), 


i«l(-h s )h s 


r dqdp 

In , 2vr 


9( ~ P 2 + (p 2 ~ j^-^q-^x'/qWx/q) . 


First, let us focus on the second term: 


h>(—h s )h s 


d ^0(-p 2 + K —) (n 2 - e rp(x x ’\ 1 il>(x'/q)i/>(x/q) 


0(\k\ - \pq\) ^ p 2 - e w{x x>) q l ^(x'/q)^(x/q) 


i+ ^ 

' c ^m-\ 

1 + 2vr q 3 

r°° r w dqdp 

Jo J-\k\ 2vr q 3 

"°° r\k\ dqd p 

/ A- 

) j-\k\ 


9(\k\ — |P|) (P 2 — k 2 ) e l i xd q /q)ip(x/q) 


— (P 2 — k 2 ) e 1 1 ( x x ^q 1 i/j(x'/ q)ip(x/q) 


j^{x/q) ( - k 2 J e P « (x x ' ] q 1 ip(x'/q) 


2 sin ^(x ~ x ')) 


~^( x /q) - q 2 ^ - k 2 --- L ${x'/q) 


2 d q / k 


n J 0 q \ q ( x - x' 


cos ( — (x — x') ) sin(^(x — x 


(x — x 1 ) 3 


ip(x' / q)'il>(x / q) 


Considering now the first term, 
4 f d( l d P ,, 


Ah = 


/ 02 k2 \ 
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dq5(x — x)q ip(x'/q) 


^,xx(x/q) + 2 ij>, x (x/q)d x + ^{x/q)d 2 x 


— k 2 / dq5(x — x')q~ ipfa 1 /q)if)(x/q ), 

Jo 

and introducing the integrals 

poo poo poo 

In-= / y n \'4>\ 2 dy , J n := / y n \ip'\ 2 dy , /C n := / y n (^' - fi^dy , 


we get the expression 


i/C 


A, = [ 1 -iP 2 + ~y(Q~ 1p + PQ- 1 ) + (C 7 i - fc 2 W 2 

This reduces to (37) if -0 is real, since then /Co = 0 and we evaluate, with the notation 

= AC-H/ACA), 

= J C 1 + ^oiH), ^ = 6i(^) 2 

The lower symbol of A s is constructed out of the following partial results: 

1 


(q,p\P \q,p) =p + t&oH&jH iy £ 32 (^0 - i)g 


-2 


and 


(g,p|Q 2 |g,p) = 60 H 62 {v)q 2 

Appendix B. Quantization of 6 

The kernel of Ag is (dimensionless) 

(64) 

« 1 


(. xlAglx') = 5(x—x')— , . 

N 1 1 ' K ’ 2ttK 0 (u) x — x' 

The semiclassical expression is 
(65) 

(q,p\A e \q,p) = 


Im 


Aik x — x' ( v x + x' 
1 - r - -K 1 - ,- 


Aik x — x' 


x + x' 


T- 


'XX' 


X + x' 


e 


dy 


ttAo(z/) 2 7 0 ^(1 - y) 


Im 


7(2/) k)Ki ( (v,k) 


Re 


An - 


i/l + i/ 


2 


7 (?/, 9P) 


with 

( 66 ) 


7(2/, a) = \ /1 


4ia 1 — y 


^ 1 + 2 / 
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